Model for a dune field with exposed water table 
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'Abstract 



Aeolian transport in coastal areas can be significantly affected by the presence of an exposed water table. In some dune 
fields, such as in Lengois Maranhenses, northeastern Brazil, the water table oscillates in response to seasonal changes 
of rainfall and rates of evapotranspiration, rising above the ground during the wet season and sinking below in the dry 
^^.period. A quantitative understanding of dune mobility in an environment with varying groundwater level is essential 
for coastal management as well as for the study of long-term evolution of many dune fields. Here we apply a model 
^ for aeolian dunes to study the genesis of coastal dune fields in presence of an oscillating water table. We find that 
the morphology of the field depends on the time cycle, T^, of the water table and the maximum height, i/w, of its 
I— I oscillation. Our calculations show that long chains of barchanoids alternating with interdune ponds such as found at 
Lengois Maranhenses arise when is of the order of the dune turnover time, whereas dictates the growth rate 
II of dune height with distance downwind. We reproduce quantitatively the morphology and size of dunes at Lengois 
O Maranhenses, as well as the total relative area between dunes. 
^ . 
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1. Introduction 



It is widely accepted that dune morphology depends 
Onfundamentally on the wind directio nality and on the 
'amou nt of sand available for transport (jWasson and Hvdd . 
^ [19831) . While longitudinal and star dunes form under 
^ bi- and multidirectional wind regimes, respectively, the 
best understood types of dune are formed by unidi- 



Mountnev and Russelll . 20091) . The dynamics of the 
water table in some aeolian sand systems is con- 
nected with seasonal variations in climate and rainfall 



( Levin et al. 



20091) ■ and has been modeled by many differ 



tectional wind. 



o 



In this case, crescent-shaped barchans 
^— ^ occur if the sand availability is low, while transverse 
'dunes, which display nearly invariant profile in the di- 
rection orthogonal to the wind, form if the ground is 
'covered with sand. The shape of dunes can be sig- 
I nificantly modif ied due to natural agents such as veg- 
'etation growt h ^Ees?, '2002'; ' Tsoar and BlumbereJ . 12002 : 
• I— I iBarbosa and D ominauez. 200 4) o r cementation of sand by 
"mineral salts ( Schatz et al.l . l2006l ). Modeling has brought 



^ many insights on dune formation in environments with 
stabi li zing vegetation fflishimori aud Tanaka, 2001; B aas , 
2OO2I: iDuran and Herrmannl. l2006al: (Yizhaa et al.l . 120071 : 



ent authors ( de Castro Ochoa and Munoz Reinosol . 1997 : 
Kocurek et aU 200 ll) . However, quantitatively little is 
known about dune field morphodynamics in presence of 
a varying groundwater level. 

One example of dunes evolving in presence of a dynamic 
water table is the coastal dune field known as "Lengois 
Maranhenses" , located in the State of Maranhao, north- 
eastern Brazil (Fig. [IJ. Long chains of laterally linked 
barchans (the so-called "barchanoids" ) extending over sev- 
eral kilometers constitute the dominant dune morphology 
at Lengois (Fig. [T]i). Dunes at Lengois alternate with 
freshwater lagoons that form in the rainy seasons and, at 
some places, nearly disappear in the su mmer. In spite of 



Nield and Baasl . l2008ai bl: iLuna et al.l . l201l[) . The influ- 
ence of sand ind uration on the shape of dunes has been 
also i nvestigated (jHerrmann et aLl . l2008l : lRubin and Hespl . 
2003). 



the field research undertaken in the past (iGoncalves et al 
20031: IParteU et all l2006l : iLevin et all l2007t iKadau et al 



The water table can also play a fundame ntal role in sed- 
imen t transport on coasts and in d eserts (iKocure^ 



19921 IRuz and Meur-Fered. |2004 iKocurek et al 



i et al 



2007 



20091 ) ■ it has remained unclear how the oscillating water 
table is contributing to the long-term dynamics and to the 
shape of dunes at Lengois. 

Understanding the effect of the water table on 
dune field evolution may help predict coastal dune 



Email addresses: marcocesarlunaSgmail . com (Marco C. M. de 
M. Luna) 

Preprint submitted to Geomorphology 



deserts ( Davidson- Arnott and Evskir*, 


19881 Enzel 


3t al. 


19991 Mountnev and Thomoson. 2002; 


Chen et all 


2004 


Grotzineer et al. . 20051 iBourke and Wrav 


, 201 1|). A phys- 



ically based model that accounts for a mathematical de- 
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Figure 1: (a) Coastal dune field named "Lengois Maranhenses" in the State of Maranhao, northeastern Brazil, located 
near 2°35'S, 42°49'W (Landsat image of May 23th of 2000). North is to the top. The landscape of the field is dominated 
by chains of barchanoids, alternating with water ponds, as those in areas I and II. (b) The arrow in the map indicates 
the location of the field. 



scription of sand transport and for the evolution of the 
terrain in response to aeolian drag and variations in the 
groundwater level is required. 

In the present work, we adapt a recently developed con- 



tinuum model for aeolian dunes ( Sauermann et al. . 2001 



Krov et al.l . 120021) to study the formation of coastal dune 
fields in an area of exposed groundwater. This model 
has been applied to study different types of dune encoun- 
tered in nature, and has proven to reproduce the shape 



surements (ISauermann et al.l. 2003c Parteli et al.. 


2006; 


Duran and Herrmannl. 


2006allbl: Parteh and 


Herr 


mann. 


20071 iDuran et al. 2008 


; Parteh et al..l2009l 


2011) 


. The 



dune model encodes a mathematical description of the tur- 
bulent wind field over the terrain with a continuum model 
for saltation — which consists of grains hopping in ballis- 
tic trajectories close to the ground and ejecting ne w par- 
ticles upon collision with the bed ( BagnoldL Il94l[) . Here 



we add to the dune model a water table that can rise and 
sink seasonally, thus affecting local rates of sand transport 
during the evolution of the dune field. In particular, we 
aim to understand the genesis and evolution of dunes in 
presence of a dynamic water level in order to shed light on 
the conditions leading to the dune morphologies observed 
at Lengois Maranhenses. 

In the next Section, we present a brief introduction 
about climate and wind regime of the area of Lengois 
Maranhenses. In Section [3] we present a description of 
the dune model and its extension in order to account for 
the water table. The calculations of coastal dune fields are 
described in Section |4l while in Section [5] we present and 
discuss our results. Conclusions are presented in Section 

E 



2. Lengois Maranhenses 

The National Park of Lengois Maranhenses is located 
on the coastal area of the Maranhao State, northeastern 
Brazil. The area comprises 155 thousand hectares delim- 
ited by the coordinates S 02°19' and 02°45', and W 42°44' 
and 43°29' (Figs. la,b). Climate at Lengois is semi-humid 
tropical with sparse vegetation, air humid ity of 68% and 
annual average temp erature about 28.5°C (jiB AM Al . 12003 : 



Floriani et al 



200 



The landscape of Lengois is characterized by the pres- 
ence of long chains of barchanoids (Fig. [5]) giving the im- 
pression of a crumpled sheet — hence the origin of the 
name: "Lengois" means "sheets" in portugese. These 
dunes extend along 75 km and penetrate inland to dis- 
tances larger than 20 km. They are detached from the 
coast by a deflation plane of width between 600 m and 
2000 m, and are migrating on top of dunes of older gen- 



erations ( Tsoar et al. . 20091) . Insights on the formative 



process and early development stages of dunes at Lengois 
Maranhenses hav e been gained from recent field research 
(jConcalves et al.l . l2003l : jParteli et all |2006[) . Sand de- 
posited by tides on the beach is eroded by the wind thus 
forming small barch ans with height betw een 50 cm and 1 m 



close to the beach ( Parteli et al. . 20061) . As the barchans 



advance downwind, they become larger, link laterally and 
give place to barchanoidal chains that can reach heights 
of 30 m. The sand of the dunes is composed by quartz 
grains of mean diameter varying between 120/im and 
350/zm (|IB AM Al . 120031 ). which are values around the aver- 
age diameter d = 250/im previously reported for grains of 
sand dunes in other fields (Ba gnold ll94ll : |Pve and Tsoar . 
19901) . 

Wind regime at Lengois is strongly unidirectional. The 
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Figure 2: Aerial photo of barchanoids and interdune ponds 
at Lengois Maranhenses. 



corresponding sand rose shows that the s trongest winds , 
which can reach velocities over 8 ni/s (| Jimenez et al. , 
1999I) . blow mainly from the East (c.f. Fig. EJi)- The 



vector length of each direction of the sand rose gives the 
potential rate o f sand t ransp ort from that direction, or the 
drift potential (|Tsoaii l200l| ). 



DP 



\u - 



Ut]f, 



(1) 



where k « 7.3, u is the wind velocity (in m/s) at a height 
of 10 m, ut « 6.19 m/s is the threshold wind veloc ity and 
/ is the fraction of time the wind was above ut (jTsoaii . 
2001r) . The resultant drift potential [RDP] is obtained by 
calculating the vector sum of the drift potential for each 
one of the directions of the sand rose. Wind directionality 
can be then quantified in terms of the ratio /3 = RDP/ DP, 
where DP is the sum of the magnitude of the drift potential 
for all directions. A value of (3 close to unity or close to 
zero means unidirectional or multidirectional wind regime, 
respectively. At Lengois, /3 is approximately 0.97, which is 
consistent with the unidirectional wind regime of the area 
( Tsoar et al.ll2009l) . 

Wind power at Lengois is negatively correlated with 
rainfall, as depicted in Fig.[3}3. During the dry season from 
August to December there is almost no rainfall and month- 
averaged values of wind velocity are much higher than 
during the wet season. From Janu ary to July, when a l 



most 93% of the rainfall takes place ( Jimenez et al. . 1999I ). 



the interdune areas are inundated due to the high water 
table thus forming ponds (Fig. [2]). The lagoons, placed 
amidst very clean sand, have no inlet or outlet and are 
exclusively filled by rain water. Their bottom is covered 
by a soft brown or gre en sheet of algae and cyanobacte- 
ria ( Kadau et al.l . l200"9h . The interdune ponds cover 41% 
of the area of Lengois. The average interdune pond area 
is 7 ha, and most of the ponds are longer in the cross- 
wind direction than in the along-wind direction. Using 
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Figure 3: (a) Annual sand rose of Aranaii (for 2003), on 
the coast of western Ceara, located at 2°49'S, 40°13'W. 
The sand rose shows that most sand-moving winds come 
from the east. The continuous lines indicate the yearly 
total wind power {DP) from the corresponding directions, 
while the dashed-dotted arrow shows the downwind di- 
rection of the resultant drift potential (RDP). The high 
value of RDP/ DP is consistent with a strongly unimodal 
wind regime, (b) Rainfall and wind power, monthly aver- 
aged for the years 1991-2007. The data were tak en from 
the NCEP/NCAR Reanalysis ()Tsoar et al.Ll2009l) . Wind 



power and rainfaU data are for 2.5°S, 40°W, and 2.9°S, 
39.4°W, respectively. 



the definition of Robinson and FriedmanI ()2002l ). the ponds 
at Lengois have a mean circularity — whe re the circu - 
larity of a sphere is one — of about 0.47 ([Levin et al 



A large fraction of the interdune lakes formed in 
the rainy season disappear in the dry season. The land- 
scape of Lengois appears to change continuously, indeed 
dune mobility implies that the lakes often reappear in dif- 
ferent places with different contours. The effect of the 
fluctuating water level on dune morphology at Lengois is 
still poorly understood. 

3. The model 

The model used in the calculations of the present work 
consists of a set of mathematical equations that describe 
the average surface shear stress (t(x,2/)) over the topog- 
raphy, the mass fiux {q{x,y)) of saltating particles and 
the time evol ution of the surface resulting from parti- 
cle transport ([Sauerniann et al.l . 2001; Krov et al. . 2002 



Duran and Herrmannl . 2006b : Duran et al. . 120101 ). Here, 
the model is extended in order to account for the water 
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table. In the model the following steps are solved in an 
iterative manner. 

(i) Given an initial topography, e.g. a terrain with dunes 
or smooth sand hills, the average shear stress field over the 
surface is cal culated solvi n g a se t of analytical equations 



developed by IWeng et al 



(|l991l ). The wind model com- 
putes, first, the Fourier-transformed components of the to- 
pographically induced perturbations in the average shear 
stress: 



(2) 




\k\ u^i) 



(3) 



where x and y are the components parallel, respectively, 
perpendicular to the wind direction, k is the wave vector, 
and kx and ky its coordinates in Fourier space; hg is the 
Fourier transform of the height profile, which is defined as 
the envelope comprising the sand landscape and the level 
of exposed water; U is the normalized vertical velocity 
profile, I is the inner layer depth of the fiow and L is 1/4 
the mean wavelength of the Fourier representation of the 
height profile; cr = y^iLk^zo/l, Kq and Ki are modified 
Bessel functions, and zq is the aerodynamic roughness of 
the surface. The average shear stress is obtained, then, 
from the equation. 



7^= To{tq/tq +f). 



(4) 



where tq is the undisturbed shear stress over the flat 
ground and tq = |fo|. 

The wind model bv lWeng et al. ( 1991 ) is valid only for 
smooth surfaces, and thus the calculation of the wind flow 
over dunes must be adapted in order to accou nt for flow 
separation at the dune brink (jKrov et al.L 120021) . For each 



longitudinal slice of the dune, a separation streamline, 
s{x, y), is introduced at the dune lee, as depicted in Fig.|4l 
The wind model is solved, then, for the envelope. 



hs{x, y) = max{/i(a;, y), s{x, y), h^}, 



(5) 



comprising the dune surface, h{x,y), the separation 
streamlines at the dune lee, s{x,y), and the water 
table, ft-w. The streamlines of flow separation define 
the so-called separat ion bubble, inside which the wind 
shear is set as zero iKrov et "al] (12002 1. The shape of 
s{x,y) is approximated by a third-order polynomial, the 
coefficients of which are calculated from the continuity of 
h, s and their respective first derivatives at the brink and 
at the reattachment point downwind, which is computed 



assum ing that s{x,y) has a maximum slope (jKrov et al 
2002h . 



(ii) Next, the local sand fiux over the landscape above 
the water level is calculated. The cloud of saltating grains 
is considered as a thin fiuid-like layer moving over the sur- 
face. Once the wind speed exceeds the saltation threshold. 
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Figure 4: Longitudinal profile of a transverse dune, h{x)^ 
and the separation bubble, s{x), both rescaled by the dune 
height, H. The shear stress, t(x), and the sand fiux, q{x), 
are normalized by the upwind shear stress, tq, and the 
saturated fiux, gg, respectively. The water level h^{x) is 
rescaled by its maximum value, ifw. Wind direction is 
from left to right. 



the sand flux increases, first, exponentially due to multi- 
plicative process inherent to the splash events at grain- 
bed collisions. However, the wind strength decreases as 
more grains enter saltation, since these have to be acceler- 
ated at cost of aeolian mom entum (Nishi mura and Huntl . 
20001: lAlmeida et al.l . l200(jl . After a saturation distance, 
the wind is just strong enough to sustain transport and 
the sand flux is maximal. By using mass and momentum 
conservation, and by explicitly accounting for the flux sat- 
uration transients, the following equation is derived for 
the hei ght-integrated mass flux of sand per unit time and 



length (jSauermann et al.L 120011 ) 



V-g-(i-k1/gs)|'?1/4, 



(6) 



{2vsa/g)ps^ul^[(u^/u^tf 
1 



- 1] is the saturated 
*t ^ — 1] the characteristic 
length of flux saturation, whereas u» ~ v^VPa is the wind 
shear velocity and pa = 1.225 kg/m^ is the air density; w*t, 



where 

flux and 4 = {'^vla/g^)/[{ut,/u 



the impact threshold (|Bagnoldl . 119411) . is about 80% the 
minimal threshold velocity (M*ft ~ 0.26 m/s) required to 
initiate saltation; g = 9.81 m/s^ is gravity, and the average 
grain velocity, Vg, is computed by takin g the steady-state 



wind velocity within the saltation lay er ([Sauermann et al 



2001: Duran and Herrmann. 2006bD. whereas a 



lave 



and 



0.43 

7 = 0.2 are empirically determined parameters 



. P 

( Sauermann et al. ■ I2001I: Id uran and Herrmannl . l2006b[ ). 

(iii) The local height, h{x,y), of the sand landscape 
evolves according to the equation 



dh/dt = -W-q/pb, 



(7) 
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where pb = 1650 kg/m'^ is the bulk density of the sand. 
If the local slope exceeds 34° , the unstable surface relaxes 
through avalanches in the direction of the steepest descent. 
Avalanches are considered to be instantaneous since their 
time-scale is much smaller than the one of dune motion. 
The downslope flux of avalanches is calculated with the 
equation, 



gavai = fc[tanh(V/i) - tanh(6'dyn) 



Whl 



(8) 



where k = 0.9 ar id ^dvn = 33° is the so-called "dynamic" 



angle of repose (jPuran et all . l2010t ). The calculation of 



the flux due to avalanches followed by the update of the 
local height through solving Eq. (O is performed itera- 
tively until the local slope is below ^dyn- The calculation 
of the avalanches applies to the whole sand topography, i.e. 
including the surface below the water level. The only dif- 
ference between the transport under and above the water 
table concerns erosion due to the action of the wind. The 
surface below the water level is obviously protected from 
erosion and can only evolve in time due to avalanches fol- 
lowing the deposition of sand incoming from dunes above 
the water level. At the interface with the exposed water, 
the sand flux vanishes and the incoming volume of sand 
is, then, instantaneously accreted to the local surface, at 
the first grid element where the sand surface is below the 
water level. 

(iv) Thereafter, the level of exposed water, h^(t), is up- 
dated. It is assumed that the wat e r tab le is roughly pro- 
portional to r ainfall (iTsoar et "aL . 20091 ). as observed in 
real situations ( Kocurek et al.L 'l992[) . The water table os- 
cillates in time (t) with a period Tw and maximum level 
i?w according to the equation: 



/iw(i) = i/wSin[27rt/rw], 



(9) 



which is chosen here to represent the nearly sinusoidal 
behaviour of average rainfall, as m onitored over s everal 
seasons, in areas of coastal dunes (| Jimenez et al. , Il999l : 
Tsoar et all . l2009l) (see Fig. [SJd)- 



The calculations are performed using open and periodic 
boundaries in the directions longitudinal (x) and perpen- 
dicular (y) to the wind, respectively. 

4. Simulations of coastal dune fields 

The initial surface is a sand beach, which is modeled as 
a flat transverse sand hill of height 1.5 m, width 80 m and 
Gaussian longitudinal profile, as depicted in Fig. [5] The 
transverse profile of the hill has small, random fluctuations 
of amplitude of the order of the grain diameter, d w 250 yum 



(jLuna et al.L 120111 ). The sand hill is subjected to a wind 
of constant average shear velocity u*. We take u» values 
within the range 0.3 — 0.4 m/s, which is t ypical for sand- 



moving winds at Lenqois Maranhenses ([Jimenez et al 



1999t IParteli et all 120061: iTsoar et all , \200m 
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Figure 5: Schematic diagram showing the longitudinal pro- 
file of the initial transverse sand hill used in the simulation. 
The upwind hatched area of the sand hill represents the 
nonerodible backshore, which has width Ax « 25 m and 
constant profile during the calculation. The maximum and 
mininum levels of the water table, both indicated by the 
dashed horizontal lines, are and — i?w, respectively. 



The calculation of a coastal dune field in the absence 
of water tab l e has been studied in previous modeling 
20101 : 



( Duran et al. 



Luna et al.L 12011 ). As shown from 



these works, the condition for the genesis of a coastal dune 
field is a saturated influx coming from the sea. In that case 
net deposition occurs at the windward foot of the sand 
s urface. The h ill does not evolve into a migrating dune 
(jKatsuki et al.L [2005.) ^ instead it remains fixed, increases 



in size and flattens (jPuran et al. . 2010t ). The fiat sand 



surface resulting from the initial hill is unstable and devel- 
ops undulations — the so-called "sand-wave instabilities" 



(jElbelrhiti et al.L 120051 ) — which evolve into small trans- 
verse dunes of height between 50 cm and 1 m. Once the 
transverse dunes reach th e bedrock, they become unstable 
and d ecay into barchans ( Reflet et al. . 2010L Parteli et al. . 
201 ll) . Average dune size increases with distance as small 



barchans merge and collide with larger, more slowly mi- 
grating dunes downwind. In this manner, a sand hill sub- 
jected to a saturated fiux becomes a source of sand for a 
field of barchans. Figure [5] shows a snapshot of the calcu- 
lation of dun e field genesis in the abs ence of water table 
(jPuran et al.L |2010L iLuna et al.Ll201ll) . 

The calculations including the water table are also per- 
formed using a saturated flux condition at the inlet, i.e. 
Qin — Is- Furthermore, the initial surface is adapted by 
including an upwind area of width Ax protected from ero- 
sion, c.f. Fig. [5l This area represents a simple model 



for the backshore of an accreting beach (jPinglen . 120051 ) . 
i.e. the zone immediately adjacent to the sea cliff that 
lies above the high-water line. At the backshore there is 
less potential to remove sand thus favouring accumulation 
(JTsoar, 2000). In fact, if the water table rises above the 
sand upwind of the crest of the hill, as illustrated in Fig. [SJ 
then, at a later time the sand surface is fully subdued by 
the water level as the hill flattens. In this case, input of 
sand ceases and no dune field forms. In order to avoid this 
problem and to make the model more realistic, we con- 
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Figure 6: Snapshots of the calculation of dune field genesis, with no water table, after ^1000 years. The calculation was 
performed with constant wind shear velocity u■^, ~ 0.36 m/s. Wind direction is from left to right. 



sider that the sand hill has a constant profile in the upwind 
backshore area, which has width Ax ~ 25 m. Whereas in 
reality the sand surface close to the beach and water table 
may have profiles much more coruplex t han in the model 
depicted in Fig. [5] ( Kocurek et al.L 120011 ). modeling of the 



detailed process leading to the accumulation of sand at the 
backshore and the formation of the sand beach is out of 
the scope of the present work. Indeed, using different pro- 
files for the sand beach or different values of Ax does not 
change the results presented in the next Section, provided 
there is a sufficient amount of sand above the water level 
that serves as an upwind source of sand for the dune field. 

5. Results and discussion 

Next, the genesis of dunes is studied by accounting for 
the presence of a water table that oscillates as described 
in Section [3] The seasonal rise and sink of the water ta- 
ble constrains dune dynamics. The evolution of the field 
depends on the quantities controlling the dynamics of the 
water table, namely the time period (Tw) and the ampli- 
tude (-ffw) of the osciUation (c.f. Eq. ©). 

One relevant quantity for dune field morphology is the 
time cycle, T^, of the water table relative to the mi- 
gration or turnover time, Tm, of the barchan, i.e. the 
time needed for the barchan t o cover a di stance approxi- 
mately equal to its own width (lAllenl . ll974l) . The turnover 
time of a barchan of height H i s given by the equation 



(jHersen et all . 12004 iDuran et all |2010| ) 



i/2 



(10) 



where = Qs/ Ph is the bulk sand flux associated with 

the wind shear velocity w*, and the pro portionality con - 
stant, a, is approximately equal to 3 ([Parteli et al.l . f2011 ). 
In the limit where Tw is much larger than Tm, dunes are 
completely leveled off during inundation since their migra- 
tion velocity is much faster than the oscillation rate of the 
water table. Sand above water level is spread throughout 
the field forming extensive flat regions in the wet season 



whereas the wind has ample time to reshape the dunes 
again the dry season. Conversely, if ^ Tm, there is not 
sufficient time to build new dunes when the water level 
sinks below the surface. The amount of sand taken from 
the dunes during the wet season is small, thus leading to 
negligible accumulation in interdune areas. Transport be- 
tween dunes in the dry season is correspondingly small, 
such that isolated dunes separated by large fiat areas dom- 
inate the morphology of the field. 

An interesting scenario occurs when becomes of the 
order of Tm. In the wet season, the lowest portions of the 
dunes, i.e. the limbs, remain under water while the dune 
crest can migrate some distance downwind, thus leading 
to a decrease in dune height. Indeed, the amount of sand 
deposited at the interdune is large enough to construct 
small dunes in the dry season. Due to their fast relative 
migration velocity, the small bedforms emerging between 
the dunes collide with the larger dunes in their front, thus 
linking the barchans laterally and leading to the formation 
of barchanoidal chains. Therefore, the typical morphology 
observed at Lengois, with chains of barchanoids alternat- 
ing with interdune water ponds, arises when the size of 
dunes is such that their reconstitution time (Tm) is of the 
order of the time cycle of the water table (Tw)- Indeed, 
taking Q w 210 m^/year, which corresponds to = 0.36 
m/s, and = 1 year, the condition ^ r„ is fulfilled 
by dunes with heights of the order of 10 m, which is con- 
sistent with the scale of dunes at Lengois. 

In Fig. [7] we present snapshots of calculations obtained 
with = 0.36 m/s and different values of T^, namely 
1 month (a), 1 year (b) and 10 years (c). Each one of 
Figs. [7^, [71d and [7J: shows one snapshot taken at a time 
where the water level is maximum, i.e. — H^j (bottom), 
and another for the case h^j = —H^ (top), where -ffw has 
the nominal value 1 m. We see that the case which leads to 
a field of barchanoids is the one in Fig.[7}D, i.e. when r„ has 
the realistic value corresponding to a seasonal oscillation. 

The morphology of the field also depends on the max- 
imum height, reached by the water. Figure [8] shows 
snapshots of simulations performed using the nominal 
value Tw = 1 year and different values of H^: 20 cm (a). 
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Figure 7: Snapshots of the calculation of dune field genesis after ~ 1000 years, obtained with wind shear velocity 
u» = 0.36 m/s and amplitude of the water table — 1.0 m. The period of the water table oscillation, r„, is (a) 0.1 
year, (b) 1.0 year and (c) 10 years. In each figure, the images on top and on bottom show calculation snapshots when 
the water level is minimum and maximum, respectively. The blue color represents the water. Wind direction is from 
left to right. 



1 m (b) and 2 m (c). As can be seen from Fig. [51 the aver- 
age size of dunes increases with distance at a slower rate as 
-ffw becomes larger. The average dune height at a distance 
of 5 km from the beach is about 28 m when = 20 cm, 
and 12 m for when = 1 m. For = 2 m, the 
dune height is roughly constant throughout the field and 
is about 4 m. Figure [3] shows the dune height at a; « 5 km 
as a function of for different values of the wind shear 
velocity, u*. We see that and play opposite roles 
for the height of dunes, i.e. while dunes become larger for 
increasing values of shear velocity, the water table reduces 
the spatial gradient of dune height in the direction of sand 
transport. 

Dune height increases with distance more slowly in the 
presence of an oscillating water table due to cyclic phases 
of destruction and construction associated with the rise 
and sinking of the water level. When the water table rises 



above the ground, sand blown from the dunes accumulates 
into wet interdune areas, leading to flattening of the dunes. 
The sinking of the water table in the dry season exposes 
the sand deposited in the lowest areas, which serve, then, 
as internal source for dunes throughout the field. The 
maximum level, -ffw, reached by the water table within 
the wet season determines the net amount of sand blown 
from dunes into flatlands, and thus the evolution of dune 
height in the field. 

Construction and destruction of dunes due to a fiuctuat- 
ing water table have been documented from field observa- 
tions at the back-islan d dune fields on Padre Island, Texas 
( Kocurek et al.l . ll992 ). These fields constitute areas where 



the genesis and formative stages of dunes in the presence 
of a water table can be observed. Early portions of the 
fields on Padre Island are reduced to a nearly planar sur- 
face during the winter, when rainfall exceeds evaporation 
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Figure 8: Snapshots of the calculation of dune field genesis after 1000 years, obtained with wind shear velocity 
u* = 0.36 m/s and time cycle of the water table = 1.0 year. The amplitude of the water table oscillation, H^, is (a) 
20 cm, (b) 1.0 m and (c) 2.0 m. The calculation snapshot shown in each case corresponds to an instant where the water 
level is maximum. The blue color represents the water. Wind direction is from left to right. 
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Figure 9: Dune height at a distance of ~ 5 km from the 
beach as a function of the maximum water level, H^, for 
different values of the wind shear velocity, u*. 



and the water table rises. Dunes are then reconstructed 
duri ng the summer when groundwater is below the sur- 



face (|Kocurek et al.l . 119921 . As on Padre Island, the early 



landscape of the dune fields produced in the simulations 
consists of 50 cm high bedforms, which are nearly lev- 
eled in the wet season and then reform in th e dry sea- 
son. As also noted by iKocurek et al. (1992), dunes in 
the "mature" stage of the field far downwind from the 
beach have more chance of survival during wet seasons 
due to their larger size. Dune behaviour depends thus on 
the maximum water level relative to the equivalent sand 
thickne ss, S, i.e. the thickness of sand if the barchan was 



leveled ([Frvberger and Deanl . Il979l) . Since the volume of a 
barchan is given by the equation V « cW^, with c « 0.017 
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Figure 10: The plot shows the relative interdune area when 
the water level is minimum (squares) and the relative area 
of interdune water ponds when the water table reaches its 
maximum height (stars) , for different values of the oscilla- 
tion amplitude, H^. The calculations were performed with 
u* — 0.36 m/s and = 1 year, and the results shown in 
the figure refer to the area between 2.8 km and 5.6 km 
downwind from the beach. The data were averaged over 
10 — 20 snapshots at time instants between 500 and 4000 
years, when the water level was larger (stars) and smaller 
(square) than 95% of and —H^, respectively. 



( Duran et al.l . 120101 ). and since V ^ 5 ■ W'^, we obtain, 

d w cW, (11) 

which can be also written, usii ig the linear r e lation between 

as (5 « 
- 10 m 



2Ql3), 



W and the dune height, H (jPuran et al. 
c- [12i7+5]. Therefore, the smallest dunes of width 
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and height '--^ 1 m are nearly leveled in the wet season when 
iJw is around 20 — 30 cm. 

The effect of the oscillation amplitude of the water 
on the morphology of the field can be further elucidated 
through Fig.[TUl This figure shows that the total interdune 
area during the dry season decreases with the amplitude 
of the water table oscillation. The larger the larger the 
amount of sand deposited in the interdune area and the 
smaller the relative area that is free of sand in the dry sea- 
son. In contrast, the relative area of exposed water during 
maximum water level increases with H^^, provided is 
not too large. Figure [10] suggests that, for intermediate 
values of i/w smaller than 1.5 m, the relative interdune 
area can serve as a proxy for the maximum water level 
relative to dune height at a given location of a dune field. 
When the value of H„ increases beyond 1.5 m, the fraction 
of dune surface above the water level becomes so small that 
erosion of dune crest leads to flat sand sheets of thickness 
comparable to the dune height. As shown in Fig.[Bl;, these 
interdune sand sheets can fill large areas of the lagoons 
between dunes, thus leading to a decrease in the value of 
/w as the water level becomes larger. Indeed, the average 
value of /w for = 2.0 m is slightly smaller than the one 
for iJw = 1.5 m, as shown in Fig. (TU] For such large values 
of i?w J the morphology of the field deviates from the typ- 
ical landscape of barchanoidal chains observed at Len^ois 
Maranhenses (c.f. Fig. Therefore, in order to make 
a quantitative comparison of the simulation results with 
dunes at Lengois, we take smaller values of -ffw, namely 
around 1 m. 

Figure lllb . shows an image of barchanoidal dunes at 
Lengois Maranhenses, which are located within area I 
in Fig. la, downwind of a recently investigated field 
of transverse dunes of heights between 7 m and 10 m 



(jParteh et all . l2006l) . In Fig. Eb, we show a snapshot of 



a calculation obtained with — 60 cm, while u* — 0.36 
m/s and = 1 year are taken in consistence with the con- 
ditions at Lengois Maranhenses. We see that good quan- 
titative agreement is found between dunes obtained in the 
simulations and real dunes. In both Figs. [TTb and lllb we 
count 14 dune crests along the longitudinal cut indicated 
in the corresponding images. Furthermore, we calculate 
the total relative area of the terrain (/w) which is covered 
by water, and find /w = 28% and 31% for the dunes in the 
image and in the simulation, respectively. 

Further downwind in the field of Lengois Maranhenses, 
in area II of Fig. la, dunes are larger and can reach heights 
of about 15 — 20 m. In Fig. [T^ . an image of barchanoids 
within area II is shown, whereas a snapshot of a simulation 
performed with = 1 m is displayed in Fig. 112b . Again, 
the calculation outcome is in good quantitative agreement 
with the image of the real dune field. The value of /w 
calculated for the real dunes is about 34%, while for the 
dunes in the simulation we obtain = 28%. 

Therefore, the relative surface area (/w) that is cov- 
ered by water in the wet season as obtained from our 
calculations is about 30%, which is close but below the 




Figure 11: (a) Image of barchanoids at Lengois Maran- 
henses, within an area of 986 m x 1964 m; (b) barchanoids 
produced in the simulation using = 0.36 m/s, = 
100 cm and r„ = 1.0 year. 




Figure 12: (a) Image of barchanoids at Lengois Maran- 
henses, within an area of 986 m x 1964 m; (b) barchanoids 
produced in the simulation using = 0.4 m/s, — 
60 cm and = 1.0 year. 



value of 41% corresponding to the total area of Lengois 
(jLevin et all . l2007l) . However, it has to be emphasized 
that the value of /w estimated from our calculations 



corresponds only to interdune ponds between chains of 
barchanoids (c.f. Figs. [TT] and [T^ . In reality, the area 
of Lengois also hosts large permanent lagoons, indeed in 
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some portions of the field the sand is distributed in an 
irregular manner leading to interdune flats of different 
sizes, which further contributes to the formation of well- 
separated barchans and a larger value of /w 

We have also considered an oscillating wind shear veloc- 
ity that is negatively correlated with the water table, i.e. 

= — £/»sin[27rt /Tw1, with Tw = 1 year, U:, » 0.30 
m/s and t/* ~ 0.1 m/s () Jimenez et al. . 1999t Tsoar et al. . 
l2009l) . However, calculations using this model for wind ve- 
locity oscillating out of phase with the water table led to 
dune patterns similar to those obtained from simulations 
with a constant wind strength. Therefore, we conclude 
that the oscillation in the water table is the relevant fac- 
tor controlling the morphology of dunes at Lengois, rather 
than the seasonal modulation of the wind strength. 



of fact, when the water table is below the sand surface, 
groundwater seeps up by capillary action to fill the pores, 
such that a subsurface layer (the capillary fringe) becomes 
saturated with water and thus protected from erosion. By 
enhancing moisture concentration and sediment cohesion, 
the capillary fringe may alter local rates of erosion and de- 
position throughout the field and thus influence the shape 
of dunes and interdune a reas (jSchenk and Frvbergeiill988 : 
Grotzinger et al.l . 120051 ). Finally, it would be interest- 



ing to combine the present model for dune formation 
in presence of a water table wit h the model for sand 



tra nsport with vege t ation growth (Duran and HerrmannL 
6at Duran et al. . 2008 : Luna et al.l , 20 111) in order to 



l2006a 

investigate the combined effect of those natural agents on 
the genesis and evolution of coastal dune fields. 



6. Conclusions 



Acknowledgments 



We extended the dune model introduced by 
mermann et al.l (|2001l ) to study dune formation in 
presence of a seasonally varying water table. The model 
was applied to investigate the genesis of dunes from 
a sand source in an environment where groundwater 
level oscillates in response to variations of rainfall and 
evapotranspiration rates. We found that the morphology 
of the field is dictated by the quantities controlling the 
dynamics of the water table, namely the time period (T^) 
and the amplitude (i?w) of the oscillation. Chains of 
barchanoids, such as those at Lengois Maranhenses dune 
field in northeastern Brazil, are obtained when is of 
the order of the dune turnover time, while determines 
the rate at which dune height increases in wind direction. 
We could reproduce barclianoidal dunes with the same 
size as those at Lengois by using realistic values for (1 
year) and (between 60 cm and 1 m) and an average 
wind shear velocity m, — 0.36 m/s that is within the 
ranges observed at the real field. Furthermore, we found 
that the maximum water level affects the total relative 
interdune area in the dune field, /w. Our calculations 
yield an estimate for at Lengois which is close to the 
value reported from measurements. 

Several factors should be included in the future in order 
to improve the quantitative assessment of dunes at Lengois 
Maranhenses. While the calculations of the present work 
considered a constant (saturated) sand infiux, an improved 
calculation of the genesis and formation of Lengois should 
consider a sand infiux which may vary in time and also 
along the transverse profile of the beach. In some portions 
of the field, transport rates may change significantly due to 
the oc curence of permane nt lagoons or areas with quick- 



sand (iKa dau et all . l2009l) . leading to spatial changes in 



sand availability. Variations of local topography may also 
contribute to the changes in sediment budget, dune height 
and local morphology of the field. Moreover, the present 
model must be improved in the future in order to account 
for the capillary fringe of the water table, i.e. the wet layer 
of sand immediately above the water table. As a matter 
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